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Abstract 


This manuscript puts forward two new generalized families of Jarratt’s 
iterative schemes for deciding the solution of scalar and systems of non- 
linear equations. The schemes involve weight functions that are based 
on bi-variate rational approximation polynomial of degree two in both its 


numerator and denominator. The convergence study conducted on the 
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schemes, indicated that they have convergence order (CO) four in scalar 
space and retain the same number of CO in vector space. The numerical 
experiments conducted on the schemes when used to decide the solutions 
of some real-life nonlinear models show that they are good challengers of 


some well-known and robust existing iterative schemes. 


AMS subject classifications (2020): Primary 65H05; Secondary 41A25. 


Keywords: Jarratt scheme, Iterative scheme, Rational approximation poly- 


nomial, Nonlinear equation, System of nonlinear equation. 


1 Introduction 


A challenging task in the area of numerical computing is that of determin- 
ing the solution of a scalar nonlinear equation I'(x) = 0 and the system of 
nonlinear equation ['(X) = 0. This is because, plethora of real life situa- 
tions or systems, are often modeled into either scalar nonlinear or system 
of nonlinear equations. For instance, More [8] presented a compiled phys- 
ical problems that are mostly modeled in the form of system of nonlinear 
equations. Problems in chemical equilibrium were modeled into nonlinear 
equations and studied by Shacham and Kehat [15]. Systems of nonlinear 
equations were applied to modeled problems in kinematics syntheses, neuro- 
physiology, economics, and combustion engineering, which were considered 
and studied by Grosan and Abraham [3]. The nonlinear global positioning 
system problems modeled into nonlinear systems, were extensively exploited 
and solved by Yaseen, Zafar, and Alsulami [19]. More literature on phys- 
ical occurrences or phenomena that are modeled into systems of nonlinear 


equations can be found in [6, 7, 9]. 


In order to have better insight of the problems that are modeled into non- 
linear equations, the solutions of the models are usually desired to be known. 
Unfortunately, there is no general algebraic formula for solving all types of 
nonlinear equations. Consequent upon this, iterative schemes are designed to 
fill this gap. While many researchers have developed several schemes for deal- 


ing with the solutions of scalar nonlinear or system of nonlinear equations, 


Tran. J. Numer. Anal. Optim., Vol. 14, No. 2, 2024, pp 391-416 


393 Jarratt and Jarratt-variant families of iterative schemes for scalar ... 


some of the schemes suffered some setbacks such as low convergence order 
(CO) and efficiency, non-convergence, divergence from true solution, break- 
down, or inability to directly extend the schemes for solving scalar nonlinear 
equations to solving multidimensional nonlinear system and still conserving 
their CO. 


The Jarratt scheme (JS) [4] presented as 


I’(2;) + 3I"(y:) ] T(x) 
GI” (ys) — 20" (aa) | T’(au)’ 


Uj41 = VU = 0,1, 252205 (1) 


where y; = Xj — 2 Ti) ig one of the famous classical iterative scheme that is 


3 I"(axi) 
efficient, converges to solution with high CO, and can directly be extended 
to multidimensional form and still retain its CO. Many JS families have been 
put forward in literature. One of such families is that due to Behl, Kanwar, 


and Sharma [1] and presented as follows: 


[(q? — 22q1g2 — 27q3)T’ (xi) + 3(qz + 10qi92 + 5930’ (ys) T(ai) 
2 [qiT’ (xi) + 3q20" (ys)] [3(qa + 92)T’ (ys) — (G1 + 5q2)T" (xi) . 


where q; € R such that q, # qo and q, 4 —3q2. Kanwar, Kumar, and Behl 


Tip. = Uy 


[5] modified the family of the iterative scheme in (2) to deal with problems 


in multivariate case. Their scheme reads as follows: 


2 
¥j = Xi — 3M(Xi) 


Kier = Xi 5 [lal + 300N(%)) Blan + a2) (X) ~ (1+ 50)" oy 


x [(at — 22q192 — 2795) + 3(qf + 10q192 + 595)N(X;)] M(X;), 
NGI), Mea On) ee, 


where J is an identity matrix corresponding to the nonlinear system dimen- 


sion. 


In [2], another JS family was designed and analyzed using complex dy- 
namics techniques to investigate whether behind the JS in the developed 
family. There are other JS variants that are computationally better in sta- 
bility and numerical performance. After analysis, their work concluded that 


JS performed best among all the concrete members of the developed family. 
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Deriving motivation from the schemes put forward in (2) and (3), we 
propose two new families of Jarratt and Jarratt-variant schemes for solving 
nonlinear equations in this paper. The new families of schemes were derived 
based on the bi-variate rational approximation polynomial of degree two in 
both its numerator and denominator. The derived families of the schemes 
are of high CO and optimal as conjectured by Traub [18]. 

This paper is arranged as follows: Section 2 is meant for the development 
of the two families and their theoretical convergence analysis. Furthermore, 
some typical members of the families are noted in this section. Section 3 
represents the extension of the families to multivariate cases and their con- 
vergence criteria. In Section 4, the application of some identified members 
of the developed families of iterative schemes to solve some real-life models 


is presented. Finally, the concluding remarks are made in Section 5. 


2 The iterative schemes 


We begin by acknowledging the following fundamental concepts. 


Definition 1. [18, 10] Let x, be the solution of [(a) = 0 and d; = |x; — z.| 
the error at ith iteration point of an iterative scheme 7(x;). If an equation 
of the form dj, = Qd7 + O(d!t') can be derived from 7(a;) via the Taylor 
expansions of the functions ['(-) and functions derivatives I’(-) it contains, 
then d;,; is known as the iterative scheme asymptotic error equation, 7 is 


CO, and 2 is asymptotic constant. 


Definition 2. [18, 10] The efficiency of an iterative scheme ¢)(2;) is the value 
nr, where T is the sum of all different functions assessment in one iteration 


cycle. 


Definition 3. [14] The computational CO (7co-) of an iterative scheme ¢)(x;) 


is the value obtained by using the formula 


log wy 


(4) 


Neoc log is ’ 


T(ai-1) 
T(aji_-2) 
three consecutive iteration results of w(«;). 


where w, = fee 


,W2= , and x;~-2,2;-1 together with x; are last 
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The formulas in Definitions 1, 2, and 3 can easily be modified to handle 


cases in vector form. 


2.1 The Jarratt’s family (JF) 


We consider the weighted scheme presented as follows: 


(zi) 
I" (a4)’ 


B41 = 04 — Poo [I’ (24), 1’) 
where 


P29 {U'(ai), 1’ (ys) 
_ @ (oi) -"/wdP + Wa C@)-Mw))+rw)) © 
ag (D" (wi) — P’(ys)]? +I" (yi) [as (wi) — P'(ys)) +P" (wi) 


is a bi-variate rational approximation polynomial of degree two in its numer- 


ator and denominator, a;(i = 1,2,3,4) € Rare free parameters such that not 


T'(2i) 


all a; = O and y = 2; - COE To determine the contribution of the scheme 


(5), some conditions are set on the parameters so as to ensure its convergence 


when implemented to solve scalar nonlinear equations. 


Theorem 1. Suppose that [ : DC R > RF is a scalar function that is 
sufficiently differentiable in D such that I’(-) 4 0 in D. Furthermore, let 
xz, € D bea simple solution of [. If 29 is in the neighborhood of x,, then 
the sequence of solution approximations {x;},.,,(@n € D) produced by the 
family of the scheme in (5), will converge to vy with CO four subject to the 


4ai—3 3—12a1+16a2 
4 16 a 


conditions that a3 = and a4 = 


Proof. Applying the Taylor’s expansion on I'(x) and I’(a) around the solu- 


tion x, and setting x = x;, then 


4 
n=2 
and 
4 
n=2 
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holds for ¢, = AVC) n> 2 


Using the expansions in (7) and (8), the expression for y; is as follows: 


3(c3 — c3) 2(4c3 — 7cacz + 34) 


f= 2; =—d; ad? de 4 d: 
ee Tig BT 5 8 3 ‘ 
fn 4(4c4 — 10c3c3 4 2 + 5C2C4 265) 45 - O(d). 
(9) 
Now, the Taylor’s expansion for ['(y;) is 
L(y) =) 4 Taig " 83 4 37¢3 d+ 2(10c3 — 16c2c3 + Dea) i 
3 9 9 27 9 (10) 
| 4 (12c5 — 27c3c3 — + 12co9c4 — 6cs) d+ O(dS), 


and for I’(y;), follow the next equation: 


2. Ach — 83 
(ys) =1+ 5d: 4 5 2a + ( 3; + Ae) d; 


3 
& 4(4c$ — Seacat 2c3 + Se2ca) 94 
(8c3 — 20c}c3 + 9eac§ ; Lc3ca ~ 3caca ~ 4cnes) d? + O(d®). 


(11) 
The use of the expansions in (7), (8), and (11), enables the deduction of 
the expansion of the the weight function P» 2 [I'’(a;),I’(y:)] in (6) as follows: 


Poo [I (xi), I’ (ys)] me 


= (a) @ 


3 
1 [2(9+ 8aq + 16a3 + 8a3 — 8a;(2 + a3) — 8a4)c3 B 
ra i 
9 + 3(-3 4 4a, — 4az)c3 


1 [4(-27 — 49a3 — 52a} — 16a3 — 4a2(13 + 4a3) 


a d? 
27 +--+ +4+27(—3 + 4a; — 4a3)c4 
1 
+ a ((8(81 + 127a3 + 21003 + 14403 + 3203 — --- + 9(4 — 5ay + 5as)cs) di 
+ O(d?). 
(12) 


The substitution of (12) into (5) yields 
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T(@:) 
=: / B / ; t 
Uj41 =X — P22 [P (a;),T (y:)| T’(a;) 
—3+4a, — 4 
=a. — dj — di 4 ( = *) @ 
3 
(Grn + 8a2 + 16a3 + 8a} — 8a; (2 + a3) — 8a4)c3 + 3(—3 + 4ay =e) 2 
9 v 
(4(—27 — 49a3 — 52a? — 16a3 — --- + 27(—3 + 4a, — 4a3)c4 Pe 
27 ss 
(Ase + 210a2 + 14403 — --. 4+ 9(4 — 5ay 4 masles ) af + O(a®). 
81 
(13) 


The expression in (13) will reduce to equation of order four if the next 


set of equations hold: 


— 3+ 4a; — 4a3 = 0, 


(14) 
9 + 8a2 + 16a3 + 8a3 — 8a;(2 + a3) — 8a4 = 0. 
The equations in (14) are satisfied when 
4a, —3 
a3 =— —_; 
4 
3- 12a, + 16a2 (15) 
a4 = ———_—__.. 
: 16 
Consequently, by the application of (15) in (13), it reduces to 
ae (S — 4a, + 16a2)et - nat) d! 
1 (4(— 1214+ 16a? + a; (48 — 64a2) — 208a2)c3 , 
108 | _ 72(-13 + 4a; — 16a2)c2c3 — 324e9c4 + 27(—-8e2 + 5) } 
+0 (a8). 
(16) 


By Definition 1, the equation in (16) is the error equation of the scheme 


in (5) and therefore, has minimum CO four. 


Remark 1. When the conditions in (15) are substituted into the scheme 
in (5) and after some simplification, a new two-free-parameter family of the 


Jarratt scheme (JF) is presented as follows: 


~, (17) 


Liga = Xj + Pog {I’(2i), I’ (ys)] Ta, 


NS 
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where P22 [I’(«;:),T’(y:)] expression is 


P29 [P" (xi), T’(y)] 


a (a [P’ (wi) — I’ (ys)]? +I" (ys) [ar (wi) — P’(ys)) + r(wi)]) (18) 
| (<8 4 1201 — 1602) I" (ay)? + 29 — 20a; + 1602) I"(a,)I"(ys) | 


+ (—31 + 28a; — 16a2I1”(y;)”) 


The JF require one ['(-) and two I’’(-) evaluations in one cycle of iteration 
and by Definition 2, the family has £.J = 1.5874. 


2.1.1 The JF Particular cases 


Here, for some specific values of a; and az imposed on the scheme in (17), a 
particular case of the developed JF can be put forward. Some typical cases 


are provided next. 


Case 1: For a, = + and az = 0, the famous JS is rediscovered as follows: 


a pects Dac 
"ar @)— 6 wi) | Py 


Vit. = VU (19) 


Case 2: For a, = # and az = 2, a new IS (JF) is put forward as follows: 


5 31" (aj)? T(a;) 
it1 = Xj t : 2 
nn =i [5+ og | Pe 0) 
Case 3: For a, = ag = 1, a new scheme (JF) is formed as follows: 
I’ (a)? — I’ (aa)T’ (ys) + I’ (yi)? T(z) 
i41 = 2,—16 . 21 
ae ee a = 100" (ai)I’ (yi) + 190" (y;)2 | T’(ai) eh) 


Case 4: For a, = a2 = 0, a new scheme (J F3) is constructed as follows: 


161" (y;)? T(2;) 
3I" (aj)? — 181" (a) (ys) + 310" (yz)? | I’ (aa) 


ig1 = Xj | (22) 
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2.2 The Jarratt-variant family (JVF) 


Here an iterative scheme that is a variant of the Jarratt’s scheme family put 


forward in (17) is presented as follows: 


Bey = ty — Roo [I (@,), 1" (y,)] at (23) 
where 
' ' _ [IP )? + bil” (@s)T’ (ys) + beI” (2)? 
Ro 2 [r (ek (yi) = l’(y,)2 ma bsI” (2; )I" (yi) ns bal" (aj)? ) 


(24) 


is a degree two rational polynomial and b;(i = 1,2,3,4) € ¥ are free param- 
eters and not all equal zero. Next, the convergence condition for the IS (23) 


is investigated. 


Theorem 2. Under the hypothesis made on the function ['(a) in Theorem 
1, the IS in (23) will be of CO four when by = ae bs = 2(1 + b)) and 


by = +57 provided by 4 —2. 


Proof. Assume that the expansions in (7)—(11) holds. Then 


I’ (a) 


R22 [I’ (zi), (ys) I’ (a3) 


1+ bi + be 
= ( ad: 
1+ b3 + b4 


1 
aH (3 + 7b3 — bo(5 + b3 — 3b4) + 11b4 + b1(—1 + 3b3 + 7b4)) a? 
2((9 + 34b3 + 33b3 + 42b4 + 90b3b4 + 6567 — b2(15 + 763 
1 38b4 — 9b] + 6b3(5 + b4)) + b1(—7 + 9b% — 6b4 + 3303 4 (2 
+— d 


963 | + b3(—6 + 34b4)))c2 + 3(—3 — 7b3 + b1(1 — 3b3 — 7b4) 


+ b2(5 + b3 — 3ba) — 11b4)(1 + b3 + ba)c3) 


1 (—4(27 + 130b3 + 23163 + 144b3 + --- — bo(19 + 2263 + --- a 
276% \ 2362(5 + 3b4) + ++» + 6b2(5 + bg — 3b4) — 11b4)(1 + bg + b4)2c4)) 
+ O(d?). 


where 06 = 1+ 63 + bg. 


Our expectation here, is to reduce the error equation in (25) to order four. 


This suffices to solving the following set of equations: 
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by + be — b3 — b4 = 0, 
4b; — b3 + 3b4 = —3, (26) 
2b; + 3b4 = —1. 


The solution set satisfying (26) is 


_ B+b, 


be 3 ’ 


i= +b); be (27) 


The substitution of (27) in (25) produces 


Ro» (l’ (xi), I’ (ys)] ae) = d; ity Gas 3b1) cS 7 3(2 ae bi )e3) 


) at + 01d 


(28) 
When (28) is substituted in (23), we have 
2 €2((10 + 3b )c3 — 3(2 + b1)cs) A 5 
it1 =x" 4 d; + O(d?). 29 
Ti41 zt ( 3(2 + bi) it (d?) ( ) 


From Definition 1, the equation (29) indicates that the IS in (23) has CO 


four. 


Remark 2. When the parameters bz and b3 are as defined in (27), then the 
scheme in (23) will reduce to a new one-parameter JVF of schemes that are 
of CO four and presented as follows: 

31" (yi)? + 3011" (x; )T" (yi) + (5 + bi )T’ (xj)? T(2;) 


BI" (yi)? + 6(1 + b1)T" (wa) (ys) — (1 + 2b1)P" (ai)? | TP’ (xa) 
(30) 


Ti4l = Vj 


Again, the family in (30) requires three assessment of functions in an 
iteration cycle. Therefore, E.J = 1.5874. 


2.2.1 Some JV family concrete members 


Here, for specific values for b;, a particular case of (30) can be put forward. 
Consider the following cases: 

Case 1: Put b; = —% in (30). A new scheme (JVF;) of CO four is 
obtained as follows: 
BI” (ai)? — 300" (ai )I" (ys) + 90" (ys)? |] Pi) 
171" (aj)? — 421” (x) (ys) + 91" (ys)? | I’ (ai) 


(31) 


Ti. = Uy 
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Case 2: For b; = —5 in (30), a new scheme (JV F2) of CO four is deter- 
mined as follows: 


I’ (yi) I" (yi) — 5I" (ai) T(z) 
I’ (ys)? — 81" (xa)I’ (ys) + 30" (xi)? | T’(ai) ” 


V4. = VW (32) 


3 Extension of the families to n-dimensional form 


In this section, the developed JF and JVF schemes were extended to solving 
the n-dimensional system of nonlinear equations [(X) = 0, where T': 0 C 
RR” — R" describes the dimension of the nonlinear system. Next, we show 
that the JF and JVF retain their CO when extended to solve n-dimensional 
nonlinear systems of equations. We define X, = (ar ae ies 2) as 
the solution of the ['(X) = 0, and for an initial iteration starting point 
Xo = (xh) 20, ia 2)" close to X,, the equivalence of the developed 
families of schemes in n-dimensional form are as follows: 


Scheme (17) (nJF): 


Xis1 = X;- A7'B[M(X,)], (33) 
where 
BI - gall — N(X))2+ 
A= | +26.N(X;i) : B= |-16] (a, (N(X;) - N(X;)*) , (34) 
+ B3N(X;)? - N(Xi)”) 


By = 12a, — 16a2 — 3, 62 = —20a, + 16a24+9, 63 = 28a, — 16a — 31, 
and 
Scheme (23) (nJVF): 


Xin. = Xi -G'H[M(X,), (35) 
where 
3N(X;)? 7 3N(X;)? 
— (1+ 2b1)f + (5+ b1)f 
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respectively. 

We note that the (¢ + 1)th iteration error of (33) and (35) is Ej4i = 
0(E7) +0 (er), is known as the error equation and 7 is the CO. Ob- 

n 

serve that HE,’ = (zP, B®), Sealy EY”) . Next, we prove that the families of 
methods in (17) and (23) retain their CO in R” space. 
Theorem 3. Assume that [: A C R” —> &” is Frechet differentiable and 
that I’(-) # O in the convex set A contains the solution X, of [(X) = 0. 
For an initial guess Xo close to X,, the scheme in (33) will form a sequence 


{Xi}; of numerical results that converges to X, with order four. 


Proof. The Taylor’s expansion of [(X) around X; up to the fourth order is 
4 

P(X) =1(Xi) +) GPO (x(x - Xi) +0 (IX- Xi’), 87) 
k=1 ~ 


where I'*)(.) is the kth-Frechet derivative of I'(-). 
Let £; = X; — X, be error at ith iteration point. Set X = X, in (37). 


Then 
4 


r%) = [Cae gr oxy By" +0 (Weal). 8) 


k=1 
Pre-multiplying the equation in (38) by [(X;)~, we obtain 


4 
M(X) = B+) [DET Ey] +0 (II). 9) 
k=2 


By applying (38) into the first step of the iteration cycle Y; = X; — 
2M(X;), we get 


2 7 i - 
Yi-Xi=a (-« a », Cue (Xi) 'T (X) (B)" +O (ve) 

(40) 

Consequently, 

(Ye — Xi)” 

= SBP - S0(%) "(BP 

+ ah (Xi) [40 (X;) + 3P"(XG)P(X;) TI" (X)] EF + O (zl) ; 

(41) 
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(Yi — Xi) =— BP + Sx) EF +0 (III), (42) 
and ie 
(¥i- X,)'=-S BF +0 (III). (43) 


The Taylor’s expression of ['(Y;) around X; is obtained as follows: 
= 2 
MY) =)>0 Face (Yi - x] +0(I%-— XP). (44) 
kai 


A little simplification on (33), while taking note that M(X;) = I’(X;)T'(X;), 
we have 
Al OD Bey, = AP BU). (45) 


Using (38) and (44), the next expansions were obtained: 
AME 

=— 160" (Xi) Ej, + (11 — day)" (X4)?T"(X4) E? 
+ ; (4(5a, — 2)(4 + a2)0" (Xj)? + (4a, — 11)1"(X4)P""(X;)) 


4 [9(28a, — 16az — 131)I’”’(X;)? + 3(92a1 — 32a2 


Bis E} 
SE = 16 Or Gar’ Gs a SP yr 
+0 (|IEill’) 
(46) 
and 
BY(X;)E; 
=— 16" (X,)E;,+ (1 — 4a, )I" (X;)°T" (X,) E? 
+ : (4(5a1 2)(4 t az)" (X;)? + (4a, = 11)0" (Xi) (X;)) 
9 [108(5a; — 4az — 6)I'"””(X;)° + 24(23a, — 8az — 38) 7 
SL sepa i (x, PO) 4 (39a, = SI’ GM ey 
+0 (||EIlI°). 
(47) 


When (46) and (47) are substituted into the right hand side of (45), we 
get 
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Gi2o, = Beg =F" CX) = P(X POR, 
AD"(X;) Buys _ 2 |6(2a1 — Bag — 5)P"(Xi) — P(X) (Xi) 
27 ie Cr (er) (48) 


+0 (ilziI) 


The error equation in (48) implies that the IS in (33) is of CO four. 


Theorem 4. Assume that the hypothesis on the function [(X) holds as in 
Theorem 3. Then the scheme in (35) will produce sequence of numerical 


results that converge to X, with CO four. 


Proof. The proof follows the manner of the proof of Theorem ??. Conse- 
quently, its error equation in R” space is 
1 | 9(10 + 3b,)P'"(X;)* — 18(2 + by 1’ ( X40" (Xi) (Xj 
GP'(X%) Ena <2 ( iI" (Xi) S 1) (XG)E" (AGP (XG) 
on 2+ br Greig) 
+0 (IEIII°). 
(49) 


4 Numerical experiments results 


This section presents the experimentation of some concrete forms of the de- 
veloped JF and JVF on the solution of nonlinear and systems of nonlinear 
equations. The numerical results obtained were compared with the results 
of some existing competitors that are also Jarratt’s variants. To verify the 
theoretical CO of the JF, JV, nJF, and nJVF obtained in Sections 2 and 3, 


we use the computational CO (7-0c) given in Definition 3. 


4.1 Numerical experiments with scalar nonlinear models 


This subsection provides some numerical experiments conducted on the de- 


veloped schemes when applied to solve some physical problems modeled into 


nonlinear equations. 
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Model application 1 (Chemical engineering, [11]). To determine the 
concentration of a chemical mixture at time x in a mixed reactor, it can be 


obtained by using the nonlinear equation 
I, (2) = —0.75 exp(—0.052) + 1 = 0. (50) 


The nonlinear equation solution is —5.7536414490.... 
Model Application 2 (Chemical Reactor, [11, 12]). In a chemical reac- 
tor, the fractional species conversion x can be determined by the nonlinear 


equation 

8x2 (4 — x)” 
(2 — x)(6 — 3x)? 
The solution to the nonlinear equation is 0.2777575428.... 


T(x) = 


0.186 = 0. (51) 


Model Application 3 (Projectile model, [17]). The nonlinear model 


governing an electron movement between two parallel plates is 


T3(x) = 7 +x2—0.5cosxz =0. (52) 


Its solution is —0.3094661392.... 

Model Application 4 (Structural engineering, [17]). The nonlinear 
model that described stresses encountered by finite structures in underground 
is: ; 

Gay z+ Eee oo, (53) 


The nonlinear model solution is 0.4160444988.... 


Model Application 5 (Population model, [16, 11]). The population 


dynamics is described by using the differential equation 


P'(x) = 6P(x) + p, (54) 


where P(a) is time x population, 6 is population birth rate, and py repre- 


sent rate of immigration. The model solution is 
P(«) = Po exp(52) + i (exp(6x) — 1), (55) 


where Pp is initial population. A problem may arise to find to find the birth 
rate 6 at the end of first year given that Pp = 100000, = 435000, and 
1564000 is population added at the end of the year. This amounts to solving 
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the nonlinear equation 


T's(6) = P(5) = = (exp(6) — 1) — 1000000 exp(5) + 1564000 = 0. (56) 


The solution of the nonlinear equation in (56) is 0.1009979296.... 


The obtained numerical experimentation results when the developed JF 
and JVF members were applied to solve the nonlinear models above were 
presented in Tables 1 and 2. The computation results in terms of the num- 
ber of iterations (IV), residual error |I'\(z)|, the total number of arithmetic 
operations with functions Ny, and computational order of convergence Negc, 
were compared with results by some existing schemes taken from [1] (see (2)) 
using qi = 0,q2 = 0 and denoted as (BS), while from [5](see (3)), taking 
q = 50,q@ = 7 and denoted as (KS). The mpmath-PYTHON was utilized 
in designing computational program and I'(z;) < 1071°° used as program 
terminal criterion. To reduce truncation error, all computations were set to 
1000 decimal places accuracy. The expression $.Te—U represents S.T x 10~/ 
where S,T,U ER. 


Table 1: Numerical results on models 1-2 


Model Is to N Nr |P(@i+1)| teoe | to N Nr |P (xi41)| Neve 
BS 5 =©80 3.3e-l11 4.0 6 96 1.3e-323 4.0 
KS 4 84 1.5e-180 4.0 4 84 1.5e-136 4.0 
JS 4 60 2.7e-178 4.0 4. 60 7.5e-134 4.0 
JF, 4 64 84e-148 4.0 4 64 7.9e-104 4.0 
JF, 4 84 5.7e-117 4.0 5 105 5.6e-295 4.0 
Ti (x) JF; 50 4 76 1.5e-1388 4.0 | 10.0 5 95 7.le-307 = 4.0 
IVF, 4 92 1.le-169 4.0 4 92 5.8e-126 4.0 
JIVE, 4. 80 9.3e-197 4.0 4 80 9.9e-135 4.0 
BS 5 80 2.4e-101 4.0 5 80 5.le-101 4.0 
KS 5 105 3.3e-269 4.0 5 105 1.2e196 4.0 
JS 5 75 5.4e-257 = 4.0 5 75 1.le-193 4.0 
JF, 5 =©80 2.2e-187 4.0 5 680 2.le-147 4.0 
JF, 5 105 1.8e-113 4.0 5 105 2.4e-103 4.0 
T(x) JF3 O01 5 95 5.8e-139 4.0 78 4 76 7.6e-141 4.0 
IVF, 5 115 = 1.7e-241 4.0 5 115  9.5e-173 4.0 
JIVE, 5 100 2.0e-139 4.0 5 100 2.2e278 4.0 
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Table 2: Numerical results on models 3-5 


Model Is to N Nr |U (aii) teoe | 0 N Nr |U (iti)! Neoe 
BS 4 64 1.6e-187 4.0 4 64 7.5e-154 4.0 

KS 4 84 7.8e-289 4.0 4 84 1.6e-193 4.0 

JS 4 60 2.0e-286 4.0 4 60 6.0e-193 4.0 

JF, 4 64 2.le-251 4.0 4 64 1.5e-179 4.0 

JF, 4 84 4.3e-219 4.0 4 84 5.8e-162 4.0 

T3(a) JF; -0.7 4 76 2.2e-250 40 |05 4 76 3.0e-222 4.0 
IVF, 4 92 6.2e-275 4.0 4 92 1.2e-188 4.0 

IVF» 4 80 2.4e-304 4.0 4 80 5.9e-222 4.0 

BS 4 64 2.6e-192 4.0 5-80 1.9e-131 4.0 

KS 4 84 9.6e-203 4.0 4 84 3.8e-106 4.0 

JS 4 60 1.2e-202 4.0 4 60 2.6e-105 4.0 

JF, 4 64 2.5e-200 4.0 5 80 1.3e-372 4.0 

JF. 4 84 4.6e-196 4.0 5 105 2.2e-307 4.0 

Ty(z) JF3 0.0 4 76 1.2e-209 40 |)09 5 95  1.7e-388 4.0 
IVF, 4 92 4.9e-202 4.0 4 92 7.0e-103 4.0 

IVF» 4 80 2.4e-202 4.0 4 80 1.5e-178 4.0 

BS 5 80 9.7e-218 4.0 6 96 2.0e-258 4.0 

KS 5 105 4.4e-391 4.0 5 105 1.7e-149 4.0 

JS 5 75 3.7e-387 4.0 5 75 2.6e-147 4.0 

JF, 5 80 6.2e-304 4.0 6 96 1.2e-415 4.0 

JF» 5 105 1.3e-227 4.0 6 126 1.2e-259 4.0 

Ts(a) JF 15 5 95 1.4e-418 40 | 30 5 95 1.7e-131 4.0 
IVF, 5 115  4.2e-351 4.0 5 115  3.7e-124 4.0 

IVF» 5 100 1.4e-132 4.0 5 100 3.8e-272 4.0 


4.2 Numerical experiments with nonlinear multivariate 


models 


This subsection put forward the experiments performed on some specific 


forms of the JF and JV to solve some nonlinear multivariate models. 


In 


all computational experiments, 100 digits floating point arithmetic and 


IP (%a) 


lloo 


<107~?° as halting criteria were utilized in the computational pro- 


gram designed in sympy PYTHON environment. Methods performance was 


compared based on the number of iterations (V), function residual error 


||\(Xi41)||, computational order of convergence ‘eoc, efficiency E, and total 
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approximate computational cost TCC. In this case, the total computational 


cost TCC of a scheme is estimated as follows: 
TCC =To4+lc+PMMct+PMVc+PSMct+MIc+ASMc+ASVo, (57) 


where P,0o,PMMc, PMVo, PSMc, MIg, ASMc, ASVc are approximate 
computational cost of evaluations in the matrix I’, vector T, product of ma- 
trix with matrix (PMM), product of matrix with vector (PMV), product 
of scalar with matrix (PSM), matrix inverse MI, addition or subtraction 
of matrices ASM, and addition or subtraction of vectors (ASV), respec- 
tively. Here, we assumed that the cost of evaluating arithmetic operations of 
functions in matrices and vectors is equal. The costs associated with these 


operations are presented in Table 3. 


Table 3: Cost of matrices operations 


Operations Cost 
I’ m2 
T m 
PMM 2m? — m? 
PMV 2m? —m 


PSM m. 
MI au 
ASM m2 
ASV m 


In Table 4, information on the approximate number of computational 
operations evaluations in one cycle of thevarious schemes are given, while 
Table 5 provides the respective scheme’s total approximate computational 


cost TCC expressed as a function of the dimension m of the system of equa- 


tions. 
Table 4: Schemes cost of matrices operations evaluations 
Methods I” [T PMM PMV PSM MI ASM ASV 

nBS 2 I I 2 3 2 2 

nk S$ 2 1 4 2 7 2 3 2 
nJS 2 1 1 2 3 2 2 2 
nJ PF 2 1 3 2 2 2 1 2 
nJ Fy 2 1 4 2 4 2 4 2 
nJ F3 2 1 4 2 4 2 2 2 
nJV;, 2 1 4 2 6 2 4 2 
nJ F 2 1 4 2 3 2 3 2 
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Table 5: Methods order and TCC 


Methods Order TCC 
nBS 4 Wm? + 11im? +m 
nKkS 4 Bn’ + 12m? +m 
nJS 4 9m3 +10m? +m 
nJ Fy 4 3278 + 6m? +m 
nJ Fy 4 23 +10m? +m 
nd F3 4 B83 + 8m? +m 
nJIV 4 28 3 +12m? +m 
nJ F 4 38m’ + 8m? +m 


MBS MKS BJS WIFI 


0.01 @JF2 @JF3 mJV1 wsV2 
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fo) 
S 
© 


0.006 
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0 | \ HTT lle mt || Balennmn 
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oO 
So 
o 
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Figure 1: Schemes efficiency dynamics for m = 2,...,5 


In line with Definition 1, we estimate the efficiencies of the schemes by 
E= NTC, Figures 1 and 2 show the dynamics of the schemes efficiencies 
as the size m changes. From Figures 1 and 2, we observe that the schemes 
efficiencies follow the order JS > BS > JF, > JF3 > JV2 > JFh > JV, => 
KS when m = 2 to m = 5, and it retains this results for the case m = 10 
to m = 30. The numerical results in Tables 6 and 7 verify these results. 
Although the compared schemes BS and KS have good efficiencies, they 
failed to converge to solutions for most problems tested. This is an indication 
that they have low stability compared to the developed schemes. 

Model Application 6 (Chemical equilibrium [7]). The following set of 
nonlinear equations ['g(X) = O describes the chemical equilibrium system 


that involves the combination of carbon oxide (x1), oxygen (#2), Hydrogen 
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Figure 2: Schemes efficiency dynamics for m = 10,...,30 


(a3), Nitrogen (4), and a5 a factor depicting the moles number of a product 


formed for each mole of consumed propane. 


1 +2122 — 3x5 = 0, 
© + 2x xq + Loxz + pal — ps + pioxd + p7T1x3 + pot2x4 = 0, 


Qxox% + 2p5x3 — 845 + p6x%3 + p7%2%3 = 0,7 


(58) 
p9t2t4 + 27 —Aps = 0, 
Ty (xo + 1) + p10Xs + Lqxs + Psx2 + psX3 
+24 —1+ pex3 + prx2x3 + pox2x4 = 0, 
where p = 10, pro = 0.193, pg = 22522, p = 2.008H8, p, — 0.00001709, py — 
eva pig = 200008846 The model exact solution X, in the domain 


A = (1,1) x (33.5, 35.5) x (-1, 1) x (—0.8, 1.8) x (—1,1) approximated to 
20 decimal places using the initial start points 
XM = (0.6,33.6,0.6,1.5,—-0.7)” and X = (—0.2,33.0,0.9, 1.0, —0.3)” is 
X, = (0.0031..., 34.5979... ,0.0650...,0.8593...,0.0369...)7. 

Model Application 7 (Combustion problem [9]). The modeled combus- 


tion problem at a temperature of 3000°C in system of nonlinear equations 


T'7(X) = 0 was given as follows: 
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—(1.0e—15) + v2 + 24% +29 + 20419 =0, —(3.0e—5) + 23 +28 = 0, 


(5.0e—5) + 21 + v3 + 2x5 + 24g +29 + 219 = 0, 


(5.0e—5) + a4 + 247 = 0, —x? + (0.5140437e—7)z5 = 0, 
— 2x? + (0.100006932e—6)x5 =0,  —22 + (0.7816278e—15)27 = 0, 
— £143 + (0.1496236e—6)xg = 0, —21X2 + (0.6194411le—7)x9 = 0, 


— £123 + (0.2089296e—14) x19 = 0. 

The solution to the model is 

X, = (1.47e-7, 2.26e—7, 1.5le—5, 6.27e—11, 
4.20e—7, 1.01e—6, 4.99e—6, 1.48e—5, 5.37e—7, 3.60e—6)", when X) = (0.1, 
0.4, 0.2, 0.3, 0.1, 0.6, 0.7, 0.5, 0.1, 0.4)" and X(? = (0.1,0.1,0.1, 0.1, 0.1, 0.1, 
0.1, 0.1, 0.1, Qa" were used as iteration starting points. 

Model Application 8 (Economics model [3]). Consider the modeling 


in economics that can be scaled up to n-dimensions and reads as follows: 


n—j—1 


Bist S- wins Inp=—-l, 1<jg<n-l. (59) 


i=l 


n—-1 
S y= -i, 
i=1 


Here, we set n = 4 to generate a system of nonlinear equations 'g(X ) = 0 
and used X() = cos (3) and X°) = (4) as starting points in two different 
applications. The model solutions is 
X, = (1.4464...,1.9730..., 4.4194... ,0.2262...)", 


Application 9 (Boundary value problem (BVP) [13]). Let the BVP 
M+ (ry?+1=0, I(0)=0, T(1)=0. (60) 


The interval [0,1] is partitioned such that % =0 < a1 < @ < +++ <a&m-_-1 < 
Im, Li41 = Vi +h, h= 4. 

Define To T'(2) 0, Ty T(21) ses | rey = T'(tm-1),l'm = T'(am) = 
1. By discretizing the problem using finite difference such that 


Dig. — Pi-1 pe Dig. — 20; +70 5-1 


, 
= 2h ‘ h2 , 


(61) 


the problem becomes 


62 
Tyga — 20; 4+Ti-1+ 7 (Tiga —Tyi) +h? =0, 2=1,2,3,...,m—1. (62) 
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For m = 10 and 6 = 2, the problem was solved by using Xp = 
{0,0,.. ao as initial start point and applying ||['(X;)|| < 10~° as stop- 
ping condition. The problem solution is 

X, = (0.605...,0.101...,0.128...,0.144...,0.149...,0.144..., 
0.128...,0.101...,0.605...)”. 

The computation results obtains on Applications 6 to 9 are presented in 
Tables 6 and 7. 


Table 6: Numerical results on Applications 6 -7 


Model Is Xo N |r (ti+1)|] Neoc TCC EB 
nBS Failed 
nk S Failed 
nJ§ 3 7.4e-07 4.0 2015 1.00069 
nJ Fy 4 1.1e-22 4.0 4287 1.00032 
nJ Fy 4 2.6e-18 4.0 5687 1.00024 
Te6(X) nJFy XM 4 10e24 4.0 5487 1.00025 
nJV Fy 3 1.4e-09 4.0 4415 1.00031 
nJV Fy 3 9.7e-08 4.0 4115 1.00031 
nBS Failed 
nKkS Failed 
nJS 4 1.7e-17 4.0 2687 1.00052 
nJ Fy 4 4,2e-19 4.0 4287 1.00032 
nJ Fo 4 2.00e-15 4.0 5687 1.00024 
Te6(X) nJFy X® 4 3817 4.0 5487 1.00025 
nJV Fy 3 3.8¢e-07 4.0 4415 1.00031 
nJV Fo 4 4.3e-21 4.0 5487 1.00025 
nBS Failed 
nk S Failed 
nJS 31 1.8¢-25 4.0 50850 1.00002726 
nJ Fy 32 1.0e-25 4.0 90123 1.00001538 
nJ FP 33 3.4e-25 4.0 122045 1.00001135 
T7(X) nJF; X 32 5.0e25 4.0 115211 1.00001203 
nJV Fy 29 2.0¢-25 4.0 110094 1.00001259 
nJV Fo 31 1.6e-25 4.0 111610 1.00001242 


4.3 Results discussion 


From the results in Tables 1, 2, 6, and 7, the new schemes put forward solved 
all the presented models with good precision and consistent CO. The eo¢ 
obtained for all the developed schemes are in agreement with the theoretical 
order of convergence derived in Sections 2 and 3. On the other hand, the 
methods of Behl, Kanwar, and Sharma [1] (nBS’) and Kanwar, Kumar, and 
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Table 7: Numerical results on Applications 7 -9 


Model Is Xo N |r (241) |] Neoc TCC EB 
nBS Failed 
nk S Failed 
nJS 29 2.20-25 4.0 47570 1.00002914 
nJ FE 31 2.20-25 4.0 87306 1.00001588 
nJ BF» 33 1.1e-25 4.0 122045 1.00001136 
T7(X) nJ F x 31 4.6e-25 4.0 111610 1.00001242 
nJV Fy 26 4.6e-26 4.0 98705 1.00001404 
nJV Fo 29 1.4e-29 4.0 104410 1.00001328 
nBS Failed 
nk S$ Failed 
nJS 9 8.8e-10 4.0 3396 1.00041 
nJ FE 6 6.2e-21 4.0 3416 1.00040 
nJ Fy 6 4.5e-21 4.0 4568 1.00030 
T's(X) nJ F3 xi 10 1.2e-22 4.0 7293 1.00019 
nJV Fy 6 2.0e-11 4.0 4760 1.00029 
nJV Fo 9 6.1e-07 4.0 6564 1.00021 
nBS Failed 
nkS Failed 
nJS 16 2.3e-18 4.0 6037 1.00023 
nJ PF 5 1.6¢e-07 4.0 2847 1.00049 
nJ Fo, 16 2.0e-17 4.0 12181 1.00011 
Ts(X) nJ F3 5 26 3.1e-09 4.0 18963 1.00007 
nJV Fy 17 1.0e-10 4.0 13487 1.00010 
nJV Fo 9 3.9e-07 4.0 6564 1.00021 
nBS Failed 
nKkS 6 1.3e-05 4.0 45252 1.00003064 
nJS 6 1.3e-05 4.0 18036  1.00007687 
nJ Fy, 6 1.3e-05 4.0 33588  1.00004127 
nJ Fy 6 1.3e-05 4.0 44289 1.00003130 
T9(X) nJ FS x) 6 1.3e-05 4.0 43308  1.00003201 
nJV Fy 6 1.3e-05 4.0 45252 1.00003064 
nJV Fs 6 1.3e-05 4.0 43308 1.00003201 


Behl [5] (n/&'S) failed to converge in most of the problems used for compu- 
tational experiments. This is an indication that the developed families of 
iterative schemes are more stable. Furthermore, results show that the clas- 
sical JS performed better than all its variants in the two developed families 
in terms of convergence and efficiency. This coincides with the findings of 
Cordero, Segura, and Terregros [2], that among the members of JS families, 
the JS is the best. 
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Conclusion 


We have successfully developed two new families of JS and its variants to 
decide the solution of scalar nonlinear and system of nonlinear equations. 
The parameters contained in the structures of the developed families can 
assume any real values subject to some simple conditions to obtain infinitely 
many special forms that have CO four when utilized to solve problems in 
both scalar and vector space. In fact, the famous Jarrat IS (JS) is a special 
form of one of the families. The computational experiments conducted on the 
new schemes show that they are good competitors to some existing schemes 
that are also JS variants. 

For further work, the complex dynamical and chaotic behavior of the 


developed schemes can be considered. 
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